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GLOSSARY OF NOTATIONS* 


Per rzontal eddy diffusivity coefficient 
Buoyancy 

Mean buoyancy 

Nondimensional buoyancy 

Puepbulenmt buoyancy fiux Cem**2/s**3) 
Drag coefficient 

Specific heat of seawater (cal/gm-C) 
Complex mean momentum 

mux Of Complex momentum 

Peompings coefficient for inertial oscillations 
Depth increment 

Total turbulent kinetic energy (TKE) 


TKE flux from atmosphere to ocean produced by 
breaking waves 


Downward TKE flux at base of homogeneous layer 
Coriolis parameter (sec*#-1) 

Downward flux of solar radiation 

Acceleration due to gravity (m/s**2) 

Rate of energy production due to mean shear 
Mixed layer depth 

Wavenumber of interfacial disturbance 


Vertical eddy diffusion coefficient for heat 
and salinity 


Vertical eddy diffusion coefficient for momentum 


Obukhov length scale 
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m,-M, Universal constants computed and used in Garwood model 
ML Mixed layer 

MLD Mixed layer depth 

D meessure 

q Square root of twice the turbulent kinetic energy 
Heat flux through sea surface 

Qn feat flux through base of ML 

Q(B) Back radiation (cal/sec*cem**2) 

Q(C) Sensible heat flux (cal/sec*oem**2) 

Q(E) Latent heat flux (cal/sec*cem**2) 

Q(N) Net surface heat flux (cal/sec*cm**2) 


Q(S) Clear sky radiation (cal/sec*om**2) 


calm Flux Richardson number 
Ri Gradient Richardson number 
Se Rossby number 
5 Salinity 
Si: Stability function for vertical eddy fluxes of heat 
amc Salinity 
Si Seabiiaty funetion for vertical eddy fluxes of momentum 


Bol Sea surface temperature (C) 


s'w' Vertical eddy flux of salinity 


t Time 
- Temperature (C) corresponding with sea surface 
" Temperature (C) 


TKE Turbulent kinetic energy 


Tw’ Turbulent temperature flux (mwatts/cm**2) 
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meno t Tone Loot ty 

Wind speed at 10m (m/sec) 

x and y components of current velocity 

x and y components of mean current velocity 
Momentum fluxes 

Z component of current velocity 


Grid~referenced horizontal corrdinates, 
positive toward east and north, respectively 


Vertical coordinate, positive upward from sea 
surface 


Reference density for water 
Density of seawater 

Melecular diffusivity coefficient 
Expansiton coefficient for heat 
Expansion coefficient for salt 
Vertical average 

Ensemble mean 

Departure from mean 

Surface wind stress (dynes/cm**2) 
Pucci pailon time scale 

Time scale of TKE transport 


Dissipation 


* Attempts have been made to maintain original notations. 
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ioe INTRODUCTION 


A. BACKGROUND AND OBJECTIVES 

Over the past few years, the ocean research community 
and the Navy have become increasingly aware of the need and 
Measures for ocean thermal structure modeling. Not only has 
the ocean response to atmospheric forcing become of impor- 
tance because of its feedback to the atmosphere, but also 
because of the importance of ocean thermal structure in 
acoustic propagation. For these and other reasons, a grow- 
ing study has centered on the upper ocean- "that layer 
bounded by the surface which exhibits strong seasonal variae- 
tions...most notable in the local temperature profile of the 
water column." [Grabowski, et.al., 1982] The nearsurface 
structure called the mixed- or mixing layer has received 
Significant emphasis. The Glossary of Meteorology defines 
this structure "as a surface layer of virtually isothermal 
water due to wind induced turbulent motion and/or free 
convection which frequently exists above the thermocline". 
To understand the form and dynamics of the mixed layer and 
the upper ocean thermocline in general, numerous models have 
come into use. 

The objective of this research is first of all to obtain 
a description of the characteristics of the mixed layer ina 


'test block' of the ocean (ca. 37 to 39N; 125 to 127W) 
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purveyed during the first and second OPTOMA cruises (Figs. 1, 2) 
occurring during the second week of March and the first two 
weeks in August 1982, respectively. (OPTOMA- Ocean Prediction 
Through Observations, Modeling, and Analysis, -is a joint 
Harvard/NPS project ‘intended to acquire field data to 
characterize synoptic scale eddies over a domain in the 
Seeeerornia Current off Northern California, and to 'set up! an 
eddy-resolving, statistical/dynamical, limited-domain, open 
boundary numerical ocean prediction model'.) The next 
objective was to evaluate the thermal structure analysis of 
the Navy's TOPS model over the second cruise period in com- 
parison with the Mellor level-2.5 model and the Garwood 

(1976, 1977) bulk mixed layer model. Finally, the results of 
the three basically single point models were compared with the 
observations. (Note: TOPS is often considered a quasi- 
three-dimensional model, because it includes advection by 
climatological geostrophic currents and wind drift.) 

The idea that the vertical distribution of temperature, 
Peiermity, and current in the upper ocean layers generally 
appears to be governed largely by the vertical heat, momentum, 
and salt fluxes imposed by local air-sea transfer [cf. 

Camp and Elsberry, 1978] has led to one-dimensional time 
dependent models which assume horizontal homogeniety. 


(Klein, 1980] The reliance upon this assumption varies 


among models as well as regions of applicability. 
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Pee TERATURE REVIEW 

Techniques for modeling vary significantly in the degree 
to which they attempt to represent the upper ocean features 
and dynamics. The existing models may be categorized as 
belonging to one of two broad classes: 

a. Bulk Mixed-Layer Models 

bee Diffusion Models 

Fach class of models will be described with some reference 
to its components. No attempt, however, will be made to cover 
all of the models in great detail or to describe each of the 
existing models or model modifications. The only exceptions 
to this will be the Garwood (1977) and TOPS [Mellor and 
Durbin, 1975] models which are of direct interest in this 
study. The structure of the Mellor Level=2.5 model will not 
be specifically discussed, because it 1s very similar to the 
TOPS (Mellor Level-2) model. 

1. Bulk Models 

The first category of models, the bulk models, are 

often referred to as slab models. The reason for this is 
that Pollard and Millard (1970) and Halpern (197%) showed, 
with the occurrence of strong wind-forced vertical mixing, 
the mixed layer responds at the inertial frequency as a 
rigid slab. This description may be misleading, for it 
implies that T, S, and V are uniform throughout the layer. 
(Glossary of Terms) Such, however, is not the case. Instead, 


the models are based on the assumption that values of the 
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aforementioned parameters are 'quasi-uniform! within a 
'well-mixed layer'. [Klein, 1980] Provided the assumption 
holds (as it may away from fronts and eddies), the values 
can be 'lumped' into integrals. [Niiler and Kraus, 1977; 
Meeacinkeyitch et.al, 1979] That is, the 'bulk' models 
predict integral (or average mixed layer) values for all 
variables. (Garwood, personal communication) -hence, the 
perhaps ‘'better' term -'integral model' is frequently used. 
This type of modeling requires an a priori assumption 
that a mixed layer exists. ([Marchuk, et.al., 1977] To 
describe the layer requires only knowledge of the fluxes 
at the surface and the base of the mixing layer (i.e., 
surface and entrainment fluxes). The mathematics of the 
modeling is, therefore, greatly simplified from partial 
differential equations to ordinary differential equations. 
Basically, mixed layer modeling is founded on two 
fundamental hypotheses. The first is that the mixed layer is 
formed as a result of the upper ocean response to forced and 
free convection. Atmospheric forcing in the form of applied 
wind stress converts mean kinetic and potential energy to 
turbulent kinetic energy within the upper ocean. This 
turbulent motion, in turn, entrains below-layer water into the 
layer, thereby changing the layer's characteristics. Strong 
solar heating works against the forced convection to 
increase the potential energy, and form a shallow turbulent 


boundary layer. This heating in effect alters only the near 
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surface water, for invoking the Second Law of Thermodynamics, 
water cannot be unmixed. Thus, solar heating results in a 
new mixed layer in addition to the ‘old mixed layer'. 
(Fig. 3) This distinction between new and old mixed layers 
Peaeeiede discussed and utilized in the following sections. 
The second hypothesis of importance, a topic of much 
controversy, invokes the use of some form of the mechanical 
energy budget as essential in the closure scheme of the more 
recent models. The disagreement, and consequently a differ- 
ence in some models, lies in the question as to whether the 
kinetic energy of the mean flow or the turbulent kinetic 
energy plays the main role in driving the downward buoyancy 
flux at the base of the turbulent boundary layer, l1.e., in 
controlling the mixed layer dynamics and entrainment. 
Thompson (1976), Pollard (1977), and Price (1977) advocate 
that the mean kinetic energy is converted directly to poten- 
tial energy by an instability in the mean flow that causes 
the bulk Richardson number to be less than some constant of 
order one. Advocates of the turbulent kinetic energy 
budget, following the original argument of Kraus and Turner 
(1967), disagree that this will be the dominant mechanism, 
reasoning that turbulence-generated instabilities normally 
have the capacity to erode the interface to the extent that 
they preclude the existence of mean flow instabilities. 


(Garwood, 1979] 
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Differences in the views by no means end there; for 
example, there is no unanimous approach to parameterization 
of the kinetic energy budget. The gross parameterization 
approach is commonly used in closure schemes. However, 
second order closure schemes have been developed, for example, 
Launder (1975, Such closure schemes, though applied in only 
one of the bulk models - Garwood (1976, 1977)-are receiving 
increasing use. Garwood's scheme will be described ina 
later section. 

2. Summary of Bulk Models 

The structure of all simple one-dimensional bulk 
models may be described in terms of two fundamental equations: 
the vertically integrated heat equation and the vertically 


integrated turbulent kinetic energy (TKE) budget. 
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where indicated terms are: 

A: Mixed layer depth (MLD) times the local time rate of 
change of mixed layer temperature (MLT). 

B: Net heat flux through the sea surface. 

See Nec heat flux through the base of the ML. 

D: P= F8n(Q. » Od 
Vertically integrated generation or dissipation of 
TKE by buoyancy forces. 

BE: Rate of TKE production due to transfer from mean 
eurrent shear. 

a ES is the downward TKE flux at the top of the ML 
(i.e., from the atmosphere to the ocean and 
produced by wave breaking at the sea surface). 


G: E, is the downward turbulent energy flux at the base 


h 
er the ML. 
Hy and H,: Effect of penetrative solar radiation. 


I: Net turbulent layer dissipation over the ML. 
J: Net local time rate of change of TKE over the ML. 
The main differences between the various bulk models 
are in the definitions of the heat flux, Q, and the turbu- 
lent entrainment across the interface. Some of the main 
features of several bulk models are listed below. 
Kraus and Turner (KT) (1967); KT developed the first bulk 
model capable of addressing the time evolution of a mixed 
layer/seasonal thermocline system; they were the first to 


peepeoze a IKE budget to close the system. Kraus and Turner 
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assumed a steady state turbulent regime and neglected 
dissipation. They parameterized penetrative solar radiation 
with an exponential decay law and shear production in terms 
mapene frictional velocity. 

Denman (DN) (1973): DN rederived the KT model, following 
Kraus! advice that the rate of increase of potential energy 
By mixing should be a small constant fraction, m, of the 
rate at which turbulent kinetic energy is transferred down- 
ward by the wind stress. (He used m=.0012) Denman also 
included absorption of solar radiation below the mixed 
layer. He showed that for sufficiently long model runs, the 
eventual depth of the mixed layer depended as much upon the 
temperature gradient below the mixed layer as upon the 
generation of turbulent Kinetic energy. Unlike KT, DN 
included some turbulent dissipation. 

momma. et.al. (PRT) (1973): As DN and KT, PRI assumed a 
steady state regime. The rate of work done by the wind 
stress on the mean motion was equated with the rate of 
increase of the energy of the mean field. The mixed layer 
deepened rapidly to a maximum within half of an inertial 
period from the time of the instantaneous increase in the 
wind speed (which then was held constant). Forty percent 
of the incoming energy was found to directly increase the 
mean kinetic energy while the rate of change of the mean 
potential energy was three times smaller than that of the 


mean kinetic energy. This result reinforced their 
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hypothesis that the shear instability of non-steady 'wind 
drift currents! is the main mechanism in turbulence 
production. They neglected turbulent dissipation and 
penetrative solar radiation. 

Nagler (NT '75) (1975): NI recognized that instantaneous 
sharp wind changes were not realistic, and that ina 

steady regime, the surface energy fluxes from the breaking 
waves could be an order of magnitude greater than the energy 
fluxes from the wind to the currents. Thus, NI extended the 
PRT model so that at the beginning and final stages of a 
model run, TKE generation through surface fluxes was most 
important, while at intermediate times the drift current 
velocity shear was dominant. As with the PRT, NI'75 

assumed a steady state, no penetrative solar radiation, and 
no energy flux through the base of the mixed layer. 
Turbulent dissipation was parameterized as a function of 

few rrictional velocity. 

Kim (KIM) (1976): KIM proposed a model more complex than KT 
which included TKE storage as well as the same parameteri- 
Zation for penetrative solar radiation as in DN. An 
additional feature of importance was the use of a background, 
depth-dependent turbulent dissipation rate to avoid infinite 
deepening. 

DeSzoeke and Rhines (DR) (1976): They used the NI model 


with the addition of a TKE storage term (term J) to show 
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that the KT and PRT models yielded results that are 
assymptotes for a single model equation. 

Gill and Turner (GT) (1976): They simulated oceanic 
seasonal cycles, and showed that the cyclical behavior of 
the potential energy content, heat content, and the 
surface temperature could be recovered if all or some 
portion of the potential energy released during the runs 
was dissipated. They neglected terms G, Hy» Hy, and J. 
Seeeber (NI '77) (1977): This model is similar to NI '75. 


Additions include a TKE storage term (term J) and background 


dissipation identical to that employed by Kim. 


Niiler and Kraus (NK) (1977): The NK model is a conglomerate 
of several of the models described above. For example, the 


radiation term is the same as that used Pim cae kh mode | 5 
terms E and F are similar to those used previously by 

Niiler with the only difference being in the coefficients; 
and the turbulent dissipation term is similar to that used 
by Resnyansky (1975) and the shear production term of the GT 
model. The NK model does not include the storage term, J. 
Unlike previous models, it includes the downward flux of TKE 
through the mixed layer base (term G). 


meeeocod (GWD) (1977): (to be discussed in Section 3) 
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SoMa POrSTERMS USED IN BULK MODELS 





er * + 
DN * + 
ERI 

ro * 75 
KIM v 
DR 

GT * % 
ha "77 # 
NK 


GWD 


* Appears explicitly 
meePppears implicitly, i.e., in one of the other terms 


Memoditied (Zilitinkevich, ¢t.al., 1979) 


3. Garwood Model 
The Garwood (1977) bulk mixed layer model is 
composed of a closed system of seven equations: 
The entrainment flux equation: 


my, <9 weiss 
e3) a Ci aaa) a 
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The horizontal component of the turbulent kinetic energy 


equation: 
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The vertical component of the turbulent kinetic energy 
equation: 
ays 
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The mean buoyancy and complex mean momentum equations: 
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The jump conditions (relating entrainment fluxes to mean 


buoyancy, rate of deepening, and mean momentum): 


—— _ oh 
(8) ~bwl=-h)=AB at ° 

— ra dh 
(9) -cw(-h)=AC ae 


Depending on which version of the model is used, 
one needs either to provide surface momentum and buoyancy 
fluxes, mean buoyancy and momentum below the mixed layer 
along with the radiation absorption or else the observed 
variables from which these fluxes can be computed. Closure 
is achieved with the bulk buoyancy and momentum equations 
and the vertically-integrated individual turbulent kinetic 
Smerey components. 

A number of significant differences exist between 
this model and other bulk mixed layer models that may not be 
obvious from the preceding summarization or the matrix of 
terms describing which terms are included in each model. 
Some of these will be made explicit in the following 
discussion of the mechanics of, and ideas beind, the Garwood 
model. 

According to Garwood (1977), the basis for the 


initialinterface destabilization and the resulting entrainment 
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is the shear across the interface (base of the mixed layer). 
This shear is hypothesized by Garwood to be the result of 
local turbulent eddies (wind generated turbulent kinetic 
energy available for mixing is explicitly parameterized as 
dependent on h/L) with the mean flow only making a minor 
contribution to the achievement of this critical shear 
meaae. (lhe instability related to this criterion is often 
Bemeea a Kelvin-Helmholtz instability or a Benjamin (1963) 
Class C instability.) By using the familiar equation for 


the flux Richardson number: 
(10) hMe=bi “uwoU/aztvwoV/ 32) 5 and 


taking into account the jump conditions (8, 9), shear 
production can be shown to be only a fixed fraction of the 
buoyant damping in the entrainment zone. This: zone may in 
fact have a Rf> 1, and still possess enough turbulent 
kinetic energy for entrainment to proceed. Garwood's 
reasoning, in contrast to that of Pollard,et.al. (1973), 
is that the mean shear production is only a secondary 
energy source for mixing. This secondary source of energy 
can only be made available by entrainment which is 
initiated by another source. This other and more 

Meat icant source is the convergence of the flux of 


Giebulent energy from above by the term: 


(1) 9 twee + BD] 
dZ 2 ome 
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Therefore, the critical number in the Garwood model is not 
the flux Richardson number, but instead the ratio of the 


Pmevyancy flux to the convergence of turbulent energy, P: 


woo pire D 


The time, T 9 


(el 3 ) eS ees 4 


required to transport the available turbulent energy, <E>, 

to the entrainment zone is important in relation to dissipa- 
tion, and to the idea that not all available TKE contributes 
to increasing the potential energy of the system. According 
to Tennekes and Lumley (1972), the viscous dissipation of 
geophysical flows with large Reynolds numbers is proportional 
to the reciprocal of the time scale of the largest eddies. 
Garwood developed the idea that for the deeper mixed layers 
(Ro~1), the planetary rotation introduces a second time scale 
of importance: (1/f). Not only did the inclusion of this 
second time scale enable Garwood to more explicitly 
parameterize viscous dissipation in relation to the local 
Rossby number, but the very role of the time scale in the 
entrainment equation enabled Garwood's model to simulate a 


Wider range of conditions, i.e., diurnal and annual ranges 


Sis, 





of mixed layer stability. (The PRT model can simulate a wide 
range of changes, though the deepening of the mixed layer is 
in accordance only with the mean shear production.) 

Another difference over most previous models was in 
Garwood's method of using the turbulent kinetic energy 
equation. The model uses the component equations (4, 5) 

Meee the typical use of the total equation. This facilitates 
analyzing the retreat or shallowing of the turbulent boundary 
layer which depends on only the vertical component of turbu- 
lence being inadequate to transport heat, momentum, and TKE 
momeenme Carlier depth of mixing. (Fig. 3) Thus, using the 
equation in component form allows for a more explicit and 
satisfactory treatment of the mixing process. 

4. Diffusion Models 

The second category of mixed layer models, the dif- 
fusion models, has received less attention from the ocean- 
ographic community than the simpler bulk models. Yet, the 
interest that exists has grown greatly with the increasingly 
powerful computers and the development of methods for solving 
nonlinear partial differential equations. ([Marchuk, et.al 
1977] The diffusion models are both complex to interpret 
and to program. They require the numerical solution of 
pertial differential equations of temperature, salinity, and 
momentum in time and depth for closely spaced vertical 
arrays of gridpoints. Unlike bulk models, the ‘'gridpoint 
models' do not require the a priori assumption of the 


existence of a mixed layer. Mixed layer depths, are, on the 
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contrary, diagnosed from the examination of the turbulence 
flux predictions as well as the predicted temperature and 
salinity profiles. 

Mathematically, the 'mixed layer’ in a diffusion 
model's vertical profile is treated in the same manner as 
is the thermocline. This allows for a more detailed analysis 
of the vertical structure (mixed layers, thermoclines, etc.), 
the mixing processes throughout the water column, and a more 
complete simulation of the vertical fluxes of temperature, 
salinity, and momentum. The latter is a basic concern in 
the diffusion models which attempt to represent the turbulent 
mixes wu’, w'T’, and w'S'’ by algebraic or differential 
equations in terms of the gradients of mean temperature, 
salinity, and momentum. These fluxes in conjunction with 
the necessary boundary conditions for T, S and u provide a 
solution for the basic conservation equations of temperature, 
salinity, and momentum. Thus, the system of turbulence and 
thermodynamic equations are closed in order to solve for 
the vertical distributions of thermodynamic and turbulence 
maepabies. (Grabowski, 1980; Kondo, et.al., 1978] 

Diffusion models may be further categorized by 
complexity, i.e., a) no-equation models; b) one-equation or 
k-models; c) two equation or k-l-models; and d) second 
Order closure models. 

A no-equation model is based on the turbulence 


transport being parameterized by a buoyancy influenced eddy 
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viscosity coefficient. The one-equation or k-models are 
Similar to the no-equation models though they include a 
relationship between eddy diffusion and the local turbulence. 
The two-equation or k-l-models include a transport differen- 
tial equation. The second order closure models, first 
proposed by Chou [1945], are based on transport equations 
for all of the Reynolds tresses derived from the Navier 
Stokes equations. They also include turbulent heat and 
salinity fluxes and the variance derived from the latter in 
conjunction with the heat conservation equation. 

A detailed profile is sought to represent the 
oceanic vertical structure as completely as possible. For 
example, Halpern (1976) among others, noted the existence of 
Shear in the top of the water column, a characteristic 
assumed not to exist in bulk models. Subsequently, the 
necessity of using numerical models to determine a more 
realistic distribution of variables such as the velocity, 
temperature, turbulence intensity, dissipation, and the like 
was recognized [Phillips, 1977]. 

The first ocean diffusion model was implemented by 
Ekman in 1905. It was a type A model, that is, it was a 
no-equation model that neglected buoyancy and salinity. A 
Notable result from this model was the famous Ekman spiral. 
Because of its neglect of density stratification, this model 
Will not be discussed further. The major difference between 


the diffusion models, beyond the number of equations 
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employed, lies in whether they are based on observations, 
dimensional analysis, and simple stability arguments or on 
some form of the TKE energy equation. Of the five models to 
be briefly discussed below, the first two are of the former 
type, while the other three are of the latter type. 
5. summary of Diffusion Models 

Munk and Anderson (MA): 1949: The MA model is a 
no-equation model wherein the eddy diffusion coefficients 
for heat and momentum are functions of the gradient 
Richardson number. No salinity is included in this model. 

Kondo, Sasano, and Ishii (KSI): 1979: The main fea- 
meme Or the SKI model is that constant flux formulations are 
applied to regions of varying fluxes. Formulation of 
diffusion coefficients based on studies of the atmospheric 
boundary layer are cast in Monin Obukhov formalism and 
applied to the upper oceans. PdG@yadisaus lon Coerficients 
for salt are assumed to be the same as those for heat. 

Vager and Zilitinkevich (VZ): 1968: This model is 
a one-equation or k-model based on Prandtl's extension of 
the mixing length theory for homogeneous flow. Salinity is 
lgnored. 

Mellor and Yamada (MY): 1974: Mellor and Yamada 
formulated a sequence of models with varying degrees of 
complexity using second order closure. This discussion will 


deal with the level-2 model, because the TOPS/TOPS-EOTS 
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model (to be discussed in the next section) is based on the 
level-2 MY turbulence closure scheme, and because it is simpler 
than the level-2.5 medel. (The level-2.5 MY model was also 
used though its mechanics are quite similar to those of the 
level-=2 MY model. The main difference is the inclusion of 
additional complexity as, for instance, the retention of the 
TKE storage term in the level=-2.5 MY model.) (To avoid 
Mepetition, TOPS/TOPS-EOTS will not be mentioned in this 
section except to say that the 'Clancy version! employed 
includes salinity, solar radiation, and horizontal advection 
by instantaneous wind drift and climatologically averaged 
geostrophiec currents.) The level-2 MY model neglects 
horizontal advection and diffusion terms. The TKE budget is 
consequently a balance between the production and dissipation 
of TKE and the conversion to potential energy as seen in 
Table II. A similar balance exists between the turbulent 
heat flux, the temperature variance, and the individual 
stress components. Unlike the other two models described in 
Table II, the MY level-2 model neglects the turbulent transport 
een . 

hMeghereim, Kiamok, and Sukhoruks (KKS): 1976: The 
KKS model is a two equation or a k-l-model. Its TKE budget 
Parameterizes the time rates of change of TKE and the 


turbulence dissipation rate. 
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TAB iE ii 
Seiten Ore tnkM>s USED IN DIFFUSION MODELS 


(Modified: Grabowski, et.al. 1980) 


TKE Prod/ Shear Buoyancy Vertical Dissipation 
dissip Produc= Conver-~ Transport of 
rate eekesg = aleiq\ (one TKE Work Done 


PReeo wn weatcalnct Pres— 
sure Gradient 


VZ 0 
MY 0 0 0 ’ 
iS * 0 


includes term 


Memeaocs Not include term explicitly 


meee LOPS/TOPS-EOTS 

TOPS/TOPS-EOTS is the Navy's new real-time ocean 
thermal structure analysis/forecast system in use at the 
Fleet Numerical Oceanography Center (FNOC). (Operational as 
of March 1983) It was developed by the Naval Ocean 
Research and Development Activity (NORDA) as part of the Navy's 
automated environmental prediction system (AEPS). The goal 
behind TEOTS is to produce reliable, accurate, real-time 
representations of the ocean thermal structure. The motivation 
for the way the model is organized (to be discussed below) lies 


in the available sources of data and the operational nature 
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of the model. Given that the upper ocean is primarily 
atmospherically forced, and the assessment that the mixed 
layer depths have proved ‘highly predictable with a variety 
of models', such as, the Denman (1973) or Mellor and Durbin 
(MD) (1975) models, NORDA has developed an operational model 
based largely on the use of the rather well-defined atmos- 
pheric variables to drive a model to describe the relatively 
data sparse oceans. Surface wind, solar heat fluxes, and 
precipitation fields are among those necessary to drive 
moira, lEOTS. 

For vertical structure, TOPS has seventeen levels 
from the surface to 500 meters, with the majority of the 
levels placed within the first 100m to enhance the vertical 
resolution where most needed. (Appendix B) The horizontal 
grid is typical of many cperational atmospheric models, 
mmen. 328i km at GON. The input fields on this grid are 
smoothed spatially by the TOPS modified EOTS. Spatial 
averaging of the fields also smooths the temporal varability 
of the fields. Of course, horizontal features smaller than 
the spacing cannot be resolved so that the capabilities are 
limited in enclosed seas, in the vicinity of fronts and 
eddies, and near boundary currents. [Clancy and Pollak, 1982] 
This last point is important to keep in mind for the 


Simulation comparisons in Chapter IV. 
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As its name implies, TOPS/TOPS-EOTS is made up of 
two components: TOPS and TOPS-EOTS. They are coupled ina 
cyclical fashion. 

TOPS-EOTS is an improved, updated version of the 
Navy's conventional objective analysis component, EOTS 
(Expanded Ocean Thermal Structure). It supplies TOPS with 
updated (real-time) initial conditions for a twenty-four 
hour forecast. Its improvement over conventional EOTS lies 
meme fact that the TOPS forecast is fed back into TOPS-EOTS 
as a first guess field for the following day's analysis. This 
not only reduces the noise resulting from the introduction of 
data into the analysis, and makes the analysis dynamically 
and thermodynamically consistent with the atmospheric 
forcing, but it also allows the model to be run for an 
extended period of time within data-sparse regions prior to 
being updated, provided atmospheric forcing data are available. 

TOPS (Thermodynamic Ocean Prediction System), the 
forecast component, is a thermodynamic model for vertical 
structure which is forced by surface fluxes supplied By 
NOGAPS (Naval Operational Global Atmospheric Prediction 
System). It uses the MY level-2 turbulence scheme previously 
described with the added effects of advection by 
instantaneous wind drift and climatologically averaged 


geostrophic currents. 
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The prognostic equations used by the model are the 


conservation equations of temperature, salinity and momentum: 
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where: 

'D' is a damping coefficient representing the drag force by 
radiational stress at the base of the mixed layer associated 
with the propagation of internal waves downward and away 


from the wind-forced region. ‘'A' is the background 


Ly 





horizontal eddy diffusion associated with the intermittent 


breaking of internal waves. [Clancy, et.al., 1982] 


There are no horizontal pressure gradients in the 


last two equations so that the u and v represent wind drift 


The geostrophic advection is provided to the u_ 
ct 


Suerent. 
terms in the T and S equations. 


the vertical 


and v 
a 
scheme, 


Within the MY level-2 closure 
parameterized 


eddy fluxes of T, S, and momentum are 
Ctaney and) Pollak (1983)) as 


Miet@lor and Durbin (1975); 


follows: 
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turbulence length scale; q=sqrt((2)*TKE); S$ 


memetlons of the bulk Richardson number: 
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Smere 1S an implied cutoff Richardson number such 
that one can determine the three empirical constants. The 


TKE budget used to obtain q is: 


ay aV Soi. = a 
21) 2qSy, Ce + oe =! ae 2qS Ho ae i57 > 0 


Sm@emeene equation for Z% is: 


The solar radiation flux 1s parameterized by using an 
extinction profile for the most common water type. (For 
the level-2.5 model, the solar radiation flux is applied 
at the surface.) 

As mentioned previously, boundary conditions are 
mainly provided by the FNOC atmospheric models. XBT profiles 


are included in the analysis when they are available. (The 
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ocean data used in the TOPS runs will be briefly mentioned 
in later sections when the model results are discussed.) 

7. Verification Testing of Mixed Layer Models 

Numerous tests have been completed or are in the 

planning stages for verification of the forecasting 
abilities of the Garwood, TOPS, and Level=-2 Mellor or Mellor- 
Durbin (MD) models as well as for the other models pre- 
viously discussed. A review of any mixed layer model 
verification testing would not be complete without at least 
mention of MILE, the Mixed Layer Experiment, which consisted 
of an intensive examination of the upper ocean in the environs 
of Ocean Station P during a 20-day period in autumn 1977. 
Excepting that brief discussion, the mainstay of the review 
of verification testing will encompass the three models 
directly involved in this study. (Though the MD model 
is essentially the core of the TOPS model, each will be 
reviewed separately, for TOPS, a first generation operational 
model, has its own inherent verification problems.) 

feeeed Layer Experiment, 19/77; The Mixed Layer Experiment 
(MILE), 1977, is a noteworthy example wherein data were 
acquired, and compared initially against simulations produced 
using the NI'77 model and subsequently other model simulations. 
Davis, et.al. (1981) used velocity profiles, measurements not 
usually available in ML verification studies in an attempt 
to not only obtain a firm understanding of the dynamics and 


thermodynamics involved, but in so doing, to also confirm 


47 





the applicability of the bulk approach to the study domain. 
M@eyetalled to acquire such a confirmation in all cases 
except those involving strong storms wherein the upper 
oceanic shear became localized a short distance below the 
Meeunermal layer. However, by accounting for the vertical 
advection and horizontal pressure gradient within the domain, 
Davis, et.al. (1981) were able to indicate a ‘reasonable 
balance! in both the heat and momentum budgets, respectively. 
Two major difficulties encountered were a) the need to 
acquire accurate data representing changes below the mixed 
layer so as to better understand the vertical distribution 
of heat and momentum within the mixed layer, and b) the need 
to acquire better expertise in being able to separate effects 
of internal waves, energetic inertial metions and internal 
tides from responses of the upper ocean to storms. [Davis, 
et.al., 1981] Despite the numerous weaknesses in the study 
described by Davis, et.al. (1981), simulations encouragingly 
seem to agree with the observations. Their 'most favorable’ 
result (Fig. 4) was for m, = 0.39 and m, = 0.48, where m4 
and m, are tuned coefficients determining ‘limits of 
excursion’ and ‘rate of response', respectively. 

a. Review of Verification Testing of Three Models 

The Level-2 Mellor model and similar models have 

been evaluated in numerous and varied simulation experiments. 
Mellor and Durbin (1975) obtained favorable results from 


tests using five weeks of data acquired at Ocean Station Papa 
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Figure 4, MILE observations versus simulations: 
NI '77, Dashed curves: observed hour 
average T(5m) and depth where T is 0.1C 
Colderseman © C5m)) solid curves: 
modeled values of mixed layer temperature 
and depth [David, et.al., 1981] 
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(SON, 145W), as well as that data acquired nearby during 
Peegee Lhe model's ability to handle the response of the 
upper ocean to an atmospheric low pressure system, i.e., 

the observed versus the predicted SST during and after 
frontal passage, was explored by Price et.al (1978) for the 
Gulf of Mexico region. The model's ability to amply predict 
the ocean's response to stronger atmospheric forcing was 
verified in the case of the passage of hurricane Eloise by 
ie coOrtunate existence of a NOAA data buoy, EB10, located 
On the track of the hurricane. Klein (1980) used the MD 
model to simulate the variability of the mixed layer depths 
in the Mediterranean Sea. The model favorably simulated 
Beers) (Clancy, et.al., 1981] the diurnal mixed layer 
response observed during BOMEX. (The Barbados Oceanographic 
and Meteorological Experiment.) 

The Garwood model has also undergone, and is still 
undergoing, verification testing, though by a more limited 
populus. One of the most notable simulations extended over 
17 years of Ocean Station Papa data. Not only was this the 
meeet simulation of its kind in length, and in the handling 
of interannual variability in mixed layer response, but it 
was also the first simulation conducted for the period of 
the spring transition. (Comparison of the hindcasts with 
the observations lends encouragement to the ability of one- 
dimensional (bulk) models to account for a large part of the 


variance on time scales from the synoptic to greater than a 
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year. [Garwood and Adamec, 1982] Additional verification 
testing includes various studies of model performance along 
the shipping track between San Francisco and Hawaii [Steiner, 
miei, and during POLEX, a component of the North Pacific 
Peeperament. [Shook, 1980] Data assimilation tests were 
performed by Warrenfeltz (1980) and by Larsen (1981). 
Be-caucemenc MY model as at the core of TOPS, one 
might say that the TOPS model has undergone extensive testing 
and evaluation. Such, however, is not the case in that 
TOPS is a first generation operational model, and as such is 
subject to inherent problems not experienced by other one- 
dimensional or quasi-three-dimensional models. The TOPS 
verification is, for the most part, a test of the quality of 
the FNOC operational data base, the source of not only the 
initial and upper boundary conditions, but also the verifi- 
cation fields. These verification fields have necessarily 
been extracted from the objective analysis or TOPS modified 
objective analysis scheme, for observations have rarely been 
available near the model's widely spaced gridpoints. The 
distribution of observations is neither uniform nor 
constant. The model is, therefore, subject to not only a wide 
range of testing conditions, but also to a greatly fluctua- 
ting source of observations. When considering the advective 
version of TOPS, one must also consider possible 


three-dimensional processes. 
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Thus, meaningful quantitative forecast verification 
seems reasonable only in regions containing an adequate 
Supply of data. The verification fields must be obtained from 
the same scheme used to acquire initial conditions, because 
different schemes applied to the same data set are likely to 
Peeauece Varying verification fields. [Barnett, et.al., 1980] 
Mmempresent FNOC daily global intake of data consists of 
approximately 200 XBT's, 2,000 'bucket' temperatures, and 
Meemeximately 20,000 satellite SST's from TIROS-N. This 
apparent lack of sufficient subsurface definition has 
mesumlted in much of the verification effort centering on the 
Biol. 

TOPS 18 an operational model, one component of a 
resource-limited naval facility model network. As such, it 
neither has number one priority in the FNOC objectives list, 
nor can it avoid being impacted by the ever evolving compu- 
tational environment. The rate of personnel change, in 
combination with the shortage of oceanographers on the staff, 
has resulted in most of the verification testing being left 
as a milestone to be achieved in the future. 

One of the major changes which has recently 
occurred at FNOC, and which will alter the workings of TOPS, 
is the replacement of the old atmospheric primitive equation 
(PE) model with the Naval Operational Global Atmospheric 
Prediction System (NOGAPS). This change took place on 


3 August 1982, midway during the model simulation runs of 
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interest in this thesis. The change was made to improve the 
higher priority atmospheric products supplied to the fleet. 
fSueemeuesy) Or the effects of NOGAPS on TOPS results is in the 
Meo tOr thé near future. It is hoped that NOGAPS will 
provide more accurate heat fluxes for ocean prediction 
analysis. 

The analysis system has evolved from OTS (Ocean 
Thermal Structure) analysis system to EOTS (Expanded Ocean 
(hermal Structure) analysis system and finally to TEOTS or 
TOPS-EOTS (TOPS modified EOTS). The last change in the 
evolutionary sequence was implemented in full on 22 March 1983- 
that is, TEOTS replaced EOTS as a bottom boundary condition 
for the atmospheric models. Previously, EOTS was the bottom 
boundary condition of the atmospheric models so that there 
was only a one-way influence between the atmospheric and 
oceanic models, 1.e., the atmospheric model drove the 
oceanic model. With the implementation of TEOTS, the feed- 
Dack loop is complete. Thus, the bottom boundary condition 
of the atmospheric model is more influenced by the physics 
Seecene situation. 

Another change to be more fully discussed later 
was a more resilient return of the model to climatology. This 
change took place on 20 September 1982. [Pollak, personal 
communication] Though not yet documented, the change 
supposedly results in the model forecast of SST remaining 


Within approximately 3C of climatology with any deviations 
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being suppressed within approximately a two-week period. (The 
scientific basis for this criterion is unknown.) (The 
climatology receives a weight in the process such that sharp 
deviations from climatology are reduced by an increase of the 
Peeencing applied to the climatological factor.) 

Most of the published verification reports describe 
TOPS model performance prior to these aforementioned changes. 
For example, Warn Varnas, et.al. (1982) performed a 60-day 
Simulation in the TRANSPAC region of the central North Pacific 
for November and December 1976. This was done using TOPS in 
memgamnection with OTS. Comparison of SST results from both 
the advective and nonadvective versions of the model revealed 
consistent skill over both climatology and persistence. 
(Fig. 6) (As will be discussed later, results also revealed 
a heat flux bias confirmed previously by comparisons of the 
heat fluxes with those calculated by Nate Clark.) Clancy, 
@et.al., (1981) performed a short time scale synoptic veri- 
fication of TOPS/EOTS (nonadvective version on Cyber 175) 
using both pattern-of-change correlation techniques 
[Dobryshman, 1972], and the apparent forecast error 
technique. The former techniquesdemonstrated TOPS' ability 
to routinely produce real-time forecasts of large scale 
changes of SST. The latter technique indicated that, with 
adequate data coverage and limited domains of strong 
atmospheric forcing, the model could achieve a skill twice 


that of persistence. 
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Hie wbaaeseyemliications Of the coupled TOPS 
mieeae TOPS/TEOTS vice TOPS linked with EOTS. TEOTS has 
proved to be less noisy and more consistent with the 
predictions of TOPS and atmospheric forcing; this result 
memeeopeciaily true in relation to the spring transition. 
Composites of forecast verification statistics for May 1981 
Mieeeace the skill of TOPS in forecasting both MLD and SST 
change out to 7/2 hours. RMS errors indicate TOPS betters 
Meeeestence out to (but not including) 72 hours for MLD, 
pe@eeonly out to 24 hours for SST. The latter deficiency 
is believed to be the result of the possible heat bias of 
Soest NOC primitive equation model. It is yet to be 


determined whether NOGAPS can resolve this problem. 


Seeeeol!UDY AREA AND ITS CLIMATOLOGY 

The domain investigated is (Figs. 1, 2) located in the 
California Current System, just south of the Mendocino 
Escarpment (Fig. 7) in the region from approximately 37.5 to 
39.6N and 125 to 127W. This system is composed of four 


currents: the California Current, the Davidson Current, the 


California Undercurrent, and the Southern California Current. 


Of main influence in the designated domain is the California 
Current, the eastern boundary current of the large anti- 
cyclonic gyre centered near the Hawaiian Islands. The 
current, an extension of the Westwind Drift (Sverdrup, 


et.al., 1942), occurs between the North Pacific atmospheric 
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high system and the semipermanent thermal low positioned 

over central California. Geostrophic and surface Ekman 
components of the flow are consistent with the temperature 
distribution, and the atmospheric high pressure system with 
its associated wind stress. An anticyclonic eddy is the 
major feature depicted by the dynamic topography computed 
momeene CCST domain. The dynamic topography for CCSII 
clearly depicts the anticyclonic and cyclonic eddies toward 
the northeast and southwest regions of the study domain which 
dominated the region during that time period and into CCSIII. 
(Reid, et.al, 1958] (Figs. 8-10) Nelson (1977) compiled a 
climatological description of the surface wind stress and 
wind stress curl over the California Current System, by 
computing the wind stress using the bulk formula in one 
degree squares from marine weather observations, primarily 
ship reports, archived by the National Climatic Center. 
Sees. 11-13) 

Nelson encountered many of the typical problems that 
arise when dealing with multi-source data as, for example, 
uneven and dissimilar sampling in space and time. Only 
approximately 12% of the observations were actually measured. 
The resulting values provide a climatological view of the 
area. The domain of interest falls within a region of 
monthly mean maximum wind stress values for August. 
(approximately 1 dyne/cm**2) The upper ocean is undoubtedly 


influenced by such forcing. Past evidence has indicated 
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Figure 9. Dynamic topography of study domain during 
CCSI (Dyn. ecm, interval 1.0, surface 
relative to 450m) 
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relative to 450m) 
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Figure 13. Long term mean August wind stress curl 
Pretds iso7=19720iNelson, 1977) 
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that the responsiveness of the coastal regime to such 
forcing. It is speculated that elements of the coastal 
response may be advected into offshore regions (and vice 
versa) by the predominance of jets within the region. (The 
small negative value for wind stress curl indicates slight 
open ocean downwelling within the study domain itself, 
whereas the coastal regions indicate strong upwelling 
tendencies.) Another point of speculation is the direct 
forcing of the subsurface eddy field by the wind stress curl 
mechanism [Mooers, personal communication, 1983]. 
Thermodynamical atmospheric forcing fields are also 
important to mixed layer evolution. Nelson and Husby (1981) 
Compiled monthly mean heat fiux fields over the California 
Current region. As with the mean wind fields, the monthly 
mean radiative and turbulent heat fluxes were computed from 
archived surface marine weather reports ranging from 
1921-1972. The marine weather observations incorporated in 
the heat flux calculations were of a highly diverse nature. 
The distribution of the observations was nonrandom and highly 
biased toward fair weather. Also problems existed due to 
sampling errors and methods of computation. Thus, large 
errors were possible in the heat flux estimates. Errors of 
#0% in each of the components of heat exchange could 
conceivably have resulted in an error of anywhere from 


10 to 70% in the net heat flux, Q(N). 
Memeets) - 0(8) - O(E) - 0(c) 
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Despite the questionable magnitude of the long term monthly 
means, however, Nelson and Husby (1981) concluded that the 
temporal and spatial consistency of the independent surface 
heat flux estimates ‘indicates that the geographic patterns 
are realistic and significant'. The wind fields were of 
Similar or better quality. 

ieeemean fields for August of Q(S), Q(B), QCE), Q(C), 
and Q(N) (Figs. 14-18), provide insight into heat gain and 
loss within the current, and more specifically within the 
study domain. Values at the domain center (38N, 126W), are: 

Mies =216.3 Watts/M**2 

Q(B)=32.77 Watts/M**2 

Q(E)=35.92 Watts/M*#2 

Q(C)=-7.24 Watts/M**2 

Q(N)=2=154.3 Watts/M**?2 

Standard errors of the means evaluated using 

SE(Q)=SD(Q)/SQRT(N) 
where SD(Q) is the standard deviation and N is the number of 
observations (253 within the study domain region), are: 

Ore) tError=3.80 (Watts/M**2) 

OCB) :Error=1.41 (Watts/M**2) 

Ore) sError=3.28 (Watts/M**2) 

OCC) -Error=1.47 (Watts/M**2) 

MON) :Error=5.20 (Watts/M**2) 
As expected, the incoming solar incidence corrected for 
Clouds and the sea surface albedo dominates the heat 


budget. The effective back radiation and the latent heat 
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Figure 14. Long term mean August incident solar 
radiation: 1901-1972 (Nelson and 
Husby, 1983] 


68 








meee soe ool (0 eo ieee BiG VS foi 123422 fll 10 119 118 f17 Lis I! 


































oq! NORA le 
eee NATIONAL MARINE FISHERIES SERVICE 
2 ‘ : PACIFIC ENVIRONMENTAL GROUP 
weve dese deen Ey er eye ce. Ee wi MONTEREY. CALIFORNIA 
47 : ae Da. 8528-1323. 43:1 .S426. 6829.22.28 124. 9934.96 
= aaa Me oe EFFECTIVE 
se ro779.2a, 6821. BZ2t. 1924. 96 Fiva2. 1495.m5. EC-UNGIA Riven 
a ae ces ‘ele Le ian diese Sanaa BACK RADIATION 
&5 nS.S81g4 523, 9822. 423.2624. 2527.5429> 34. 1341.30 { WATTS/Mes2 3 
So be re ee ee toe es ie 
Ad ; | NG. PH28. 4424. 2521. 2623. tcc We pi 
OSE eerie pio): 0 aaa eee i a laa LONG TERM MERN 
: : ; Seeds 23. 26Z7.6521.7925. 2 P4- 3 CAPE BLeNCD FOR 
wee edwee oda sten d-=- 2 © FOOD SOE IOS TOO SOG wae dans. 
y| 2 22.6824. S626. 2225. 29. 71.20. 2525. 97 RUGUST 
~ ood Scope cee gees 2S O0D000S9 Ro acho oe odeces as ~ eee OTs. oo lee Ph 2c =i 
4] : ae PS. 1423. WS. WZ7. W24. 873A, BIA, BI 3A. 72.44, 16.39. 36- 
ae nes we wade een doccedonee sone. se. Tyran aC ae ae A aN "- 7 CAPE MENOOC [NO 
A+ ti Be. 622. 7027.8621. 7823. 2827. pa 37.7248) 3930.75. 
weeeendewwe tec oe dan = = + Goren Yorn edea edn ean da e+ de - -Fdoward-- oe tod. = 
QA + + %  p3.ge2s.2241.9325.7125. 4928. 1836.56.98. 8239.22 45. EP | 
eeeaedeeceedeowecadeoceces secet- ete- pedo rood ---e - ose & ae 
a | ae D7. 628.2426. 7728. 1728.5728. e432. 7737.52 34. 2038. EL 
SOI SRZOOB OTS GOS DOR OO HIG. eocedorce ds we- 4- - ee 4n- ww Gor e- Here eot- we: 4s -- 
37 ; of $f) 3 b.ag27. 2827.77 300zaz9. 2431.56.36. 93402535.2735.18 47 
3g > | A y 2326. 7625. 925. 623}. (32.8131.5234. 3138. 3 
POD BOSS CHO00 OBI GOO MOS OS CC. beset Sore Ome v9 eg dens 
rr 5.2324, 6924.01. of pa it ae PRE cor és 5 
rr aes fo za( ee. 5427. 19245 5728,5227. 
3 | sae ne eer a , sl 4 .3727.8221.7 6.05 73.2425. 2328. 1628: 
2 a aa a a a alas a lees (oo eG eee os ee 
3 I ee ce on ne a no “422. 7321. 7222. 3228. seas 2627. 23 
31 Re PS, 9826, 4626. 4121.64 17.56 18. 921. 3326. 2523. 7526.29 
O 2 ae ; ; 7¥2.24. 5722. 8929409 16. 1519 Gi 22.6421. 13.24.6323. 
loose: iced a OOOO : ae 
< SS 1S.372). 0528. 36:16. 27 33. 2883. 6424. S326. 9895. 
Scie Aninlels\ animal mnimaldaicinini tanned mininiai ta wnindann\odwacicdawe=(Jinccin de cca o-- 5 %-22-4--- 5 dae Oe a oe ae rae eect. Gopey boc 
vacuole re es ee Se eR &2 sar OE "SS tase acme Puna Eucent 
TS Se ee a a a PL. 2226. 2427. 623.3721. 97 2B. 4921.56.31 31 a7ab, 845.0 
1! ’ ' ‘ 1 ‘ 1 ! ! ’ ‘ 1 1 1 7 an \ 
D ' ! ' ' ' ' ’ ; ' 7. a HIB A828. 2122. 3524. 6326.21. 
= Se a ee a a 6. 4917. 26.351 (Tal .6823. 2924.81.28. 
a i ala ah a al fetta ee eee 4 
Pood bo bob a ob bog ob of ob ob bb peaebz sel zen zazsne 
PS I a aS aa ia a ae ld Sar aria ae (oleate (a See oe 
2 en a PE.87 22. NG. 29.19 S526. 5729. 14.24.5226. 7434. 4134, 
Re a a hg. 727.691 «1625. 4626, 0226. 7728. 6229. Yai. 4330. 780 
Be ala a ch i dia aa ar fl tae | al. chat Co Naa el creel 
21 ee a Ee hg. 38 TFA. Lt $7128.11 23.878. 3930, 5230, 2231. 15.2! 


peeeoms sess bog s2 1h) R129 128927 126 Az5 eles 12a te) 128 119 16 fi7 tte 11S 114 tis 


Figure 15. Long term mean August effective back 
radiation: 1901-1972 [Nelson and 
Husby, 1983] 
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Figure 16, Long term mean August sensible heat 


flux: 1901-1972 [Nelson and Husby, 
1983] 
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Figure 17. Long term mean August latent heat 


Flux: 1901-1972 [Nelson and Husby, 
1983] 
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Pigure 18. Long term mean August net heat 
exchange: 1901-1972 [Nelson 
and Husby, 1983] 


a2 





term are essentially of equal climatological value within 
the region of interest. As will be seen in the following 
@iscussions, the climatological value of net heat flux is 
far greater than that computed from acquired synoptic 
atmospheric data due to frequently overcast conditions. 

The climatological Q(N) is, however, somewhat less 
(approximately 30 Watts/m**2) than the value 'backed out' 
Meemerne transition depicted in the acquired TOPS/TOPS-EOTS 
profiles. Heat fluxes and wind stress values, the boundary 
conditions required for the model simulations, are to be 
further discussed in Chapter IV. 

Wind speed cubed, thermocline strength, and mixed layer 
depth fields are also included here to demonstrate the 
spatial and seasonal characteristics of wind generated 
turbulent kinetic energy and the associated upper ocean sta- 
bility based on Husby and Nelson (1982). Winds used in the 
computations of the wind speed cubed fields were the same as 
winds used in the calculations of wind stress and wind 
stress curl. The wind speed cubed fields were computed as an 
indicator of the turbulent kinetic energy transferred from 
P@emdir to the upper ocean, a common procedure. 

ieepmoducing the other two fields, Husby, et.al. (1982) 
were faced with some of the typical problems one encounters 
when working with mixed layers, more specifically with 
Bweerpreting temperature profiles (mainly XBT). It was 
necessary to establish the criterion for the extraction of 


the desired features. A typical vertical temperature profile 
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Seemecnaracteristics of interest is depicted in Fig. 3. (This 
Meerrle and some actual XBT profiles will be discussed in 
Seema etalil in Chapter III, Section A.) That BT profiles 
comprising their data base were digitized at five meter 
intervals plays an important role in the type of criterion 
used. They chose to select their mixed layer depths as the 
Meperstimit of the first five meter interval in which the 
Remperature gradient exceeded -0.3C/5m (-0.06C/m). The 
Dottom of the thermocline was obtained by choosing the upper 
depth of the two successive five meter intervals within 
which the temperature gradient did not exceed -0.3/5m. The 
Strength of the thermocline was then obtained by taking the 
difference of the temperature at the bottom of the thermo- 
cline from the mixed layer temperature. These values were 
then summarized into one degree squares and mapped. (Fig. 19) 
The practice of using the wind speed cubed as an indica- 
tor as described above does not seem to apply to the study 
domain in August, for in using an above average estimate of 
6 m/s wind, the resulting wind speed cubed is 216 (m/s)**3; 
a value quite a bit less than that depicted in the climatology. 
The mean observed mixed layer 1s, however, quite consistent 
with the corresponding climatological mixed layer depicted in 
Fig. 19. The discrepancy between turbulence indicator (wind) 
and mixed layer depth leads one to consider the probable 
importance of convective activity and horizontal advection 


within the region. 
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A. DATA ACQUISITION 

The analyzed data were acquired aboard the R/V ACANIA 
Semetipport Of the Pilot Ocean Prediction Study of 
California Current Eddies, which has since become known as 
the OPTOMA program sponsored by ONR. The cruises designated 
CCS1 and CCS2 (Legs I and II) covered time intervals 8 to 
Meme, 30 July to 5 August, and 9 to 14 August 1982, 
respectively. (For the convenience of this thesis, CCS1 
fees Called CCSI, CCS2 legs I and II will be called CCSII 
mmc colil, respectively.) Each had as its domain of 
interest, a region off Point Arena centered at approximately 
38N, 126W. CCSII was devoted to sampling the largest domain, 
approximately 180 km square with line spacing of ca. 44km. 
and station spacing of approximately a quarter of the 
Daroclinic Rossby radius of deformation, 1.e.,9 km. This 
allowed for better resolution over a scale neither so large 
as that of the CalCOFI data nor as small as that used during 
Meet, (120 km with line spacing of ca. 18 km and station 
meeeine of 9 km.) CCSIII was for the most part a resampling 
of the CCSII grid. Station spacing on a whole was maintained 
at small enough intervals (ca. 9 km for August cruises) to 


minimize aliasing by the smaller scale activity. 
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Meeell three cruises, continuous sampling was achieved 
Beane the shipboard DAS (Digital Data Acquisition System). 
T-4 XBT's (expendable bathythermographs) maximum depth of 
Memeeeevere launched at each station. (Figs. 1, 2) 189 XBT's 
were deployed on CCSI, while approximately twice that number, 
393, were deployed in the combined latter two cruises. No 
CTD casts were made during CCSI. Interspersed along the 
meade on CCSII and CCSIII were 500 and 1500m CTD casts. 
Meee CCSIT and CCSIII, 100 CTD casts were made prior to 
equipment failure due to rough weather. In addition to 
logging these vertical profiles, DAS stored station and 
interstation "underway" values of position, time, ship's 
Beeea, Ship's roll, and surface atmospheric data; such as, 
dewpoint temperature, air temperature, barometric pressure, 
wind direction, and wind speed. Fractional cloud cover was 
recorded in the ship's log on an hourly basis by the deck 
watchstander. 

Winds for each day of the August cruises in the vicinity 
of the cruise tracks (nine gridpoints) were supplied by FNOC. 
Meee, 20, 21) In addition, TOPS and TOPS-EOTS produced 
ocean temperature profiles associated with eight of the nine 
points of interest which were archived by FNOC for later 
Seeruation with in situ data. (The ninth point was located 
on land.) 

Satellite IR imagery was available from NOAA/NESS 


(Redwood City) for cruise planning and preliminary analysis. 
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Peecure 20. COSI Cruise track and the nearest 
TGEs/ LOPS=2O0 lomeraa points 
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-- Cruise track 
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Figure 21. CCSTIIITcruise track and the nearest TOPS/ 
TOPS-EOTS praidpoints 
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Higher quality imagery was supplied by SOSF (Scripps Ocean 


Meeeeate Facility) for post analysis. 


B. DESCRIPTION OF DATA AVAILABLE 


im, CCSI 





The XBT data are of primary interest from the March 
cruise. 189 XBT's were deployed, of which 180 were 
‘successful’. Each trace was digitized by FNOC because 
DAS had malfunctioned. Final selection left 105 XBT's 
Meee the grid region. Data files consisted of date, time, 
latitude, longitude, and vertical profiles of depth versus 
the corresponding temperature and climatologically-based 
fewanity . 

meee Coll and CCSITIT 

These cruises contained more data of present interest. 
Memo) Count was 257 after editing. Digitization by FNOC 
in this case was not necessary. Of these XBT's, 198 were 
located within the grid. 

The data, similar to that of CCSI with the exception 
that the salinity was not climatogically-based, were recorded 
on DAS cartridges. Data from these tapes were subsequently 
made accessible from CMS (Coversational Monitor System) on 
the NPS IBM 3033. The original traces and XBT logs were used 
for comparison and error checks. 

DAS log data, stored and accessed in the same way as 


the XBT data, though in the standard DAS header file format, 
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feemene Source of the majority of the atmospheric data 
memmered during CCSIIT and CCSIII. Numerous runs in various 
editing modes were made to delete blocks of zeros, erroneous 
oes, and data not applicable to the grid region. (see 
Appendix A) 

Additional atmospheric data, such as cloud cover, 
fee@eeeccorded on an hourly basis in the ship's log. Any 
additional data pertaining to the ocean, such as bucket 
Mem@peratures, were logged in the CTD log or on the XBT 
~erdaces, 

Finally, atmospheric as well as oceanic analyses in 
the vicinity of the experimental domain were archived by FNOC. 
They consisted of daily wind values, one for each of the 
Mine selected FNOC gridpoints, as well as daily analyses and 
24-hour forecasts of the ocean temperature profiles for the 
Si2nt gridpoints located over the ocean. Since TOPS/TOPS-EOTS 
runs on a Northern Hemispheric polar stereographic grid with 
@eerid spacing of approximately 300 km in the midlatitudes, 
it was not possible to obtain analyses from numerous points 
Maethe cruise region. In fact, of the eight points located 
over the ocean, none of them fell directly within the cruise 
@emain. (Figs. 20, 21) The use of the standard FNOC gridpoint 
was an unnecessary limitation that was imposed for convenience. 
(Subsequently, it was learned that FNOC has the capacity to 
interpolate their analyses to produce a vertical temperature 


profile at a specified point.) The analyses and TOPS 
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predictions used in this study were from the point closest 
Semegeestudy domain. Based on the SST climatology (Fig. 22), 
Sgemenosen FNOC gridpoint at 39.05N, 127.75W was not 
expected to have a SST much greater (at most 1C) than that 
of the XBT's within the cruise area. 
Seer oRncING FUNCTIONS FOR THE GARWOOD AND MELLOR LEVEL-2.5 
MODELS 
The version of the Garwood model utilized required an 
initial vertical temperature profile as well as three-hourly 
values of air temperature, dewpoint temperature, wind speed, 
ememmeroud Cover for the computation of the boundary conditions 
to be used in both the Garwood and the Mellor 2.5 models: 
the momentum and heat fluxes. 
mer che first run, initialization was done with the FNOC 
analyzed ocean temperature profile at 39.0N, 127.75W. This 
initialization permitted comparison of the forecasted TOPS/ 
TOPS-EOTS ocean temperature profile with that obtained from 
the Garwood model over the same time span. All other runs 
were made using XBT 96 as the initial temperature profile. 
The three-hourly air temperature was obtained by 
averaging the DAS underway data. Averages were made over 6 
to 299 values. The seeuareta standard deviations range from 
Semen Of 1.657C to a low of 0.0005C. Only three values had a 
standard deviation over 1.0C. Averaging was done under the 


assumption that the air temperature varied little over the 
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West Coast Semi-monthly means of sea surface temperature 
15 DAY SST MEAN (°C) compiled by the National Meteorological Center 
of the National Weather Service. : 


Figure 22. Average SST for last half of July 19382 
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Meee vely small cruise grid at any given time. In addition, 
it was assumed in the case of missing or erroneous data 

that, at a given time of day, the values did not vary much 
me -eecne Course of a couple of days. Thus, for missing 

eee 0300 readings, the 0300 input to the model was 

Beeumed equal to that of the previous day. (Fig. 23) 

The dewpoint temperatures were to be acquired in a 
Similar manner, however, faulty instrumentation precluded 
this. Dewpoint temperatures were instead computed using a 
Memmerguc described by Bolton (1980). (Fig. 24) To use the 
method, relative humidities (also unavailable) were required. 
As recommended by Davidson (personal communication), a 
relative humidity of 90% was assumed throughout to correspond 
with the foggy/hazy conditions which existed intermittently 
through the cruises. 

Wind speeds (Fig. 25) from the nearest FNOC point were 
used consistent with the location of initial values. The 
region of interest was considered to be too small to require 
interpolation of winds. (Fig. 26) Ship winds available from 
DAS are relative winds, but true winds could be calculated 
Dy correcting for ship's velocity. This was not necessary 
for the runs described. Instead, the applicability of the 
FNOC winds was checked by comparing them with the DAS winds 
recorded at CTD stations because the ship was then 'stationary'. 


(see Appendix A) 
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sea surface temperatures were obtained from individual 
XBT's along the track. (Fig 27) The XBT's used were those 
launched within half an hour of the necessary input time. 

m@ifesfinal input, cloud cover in octals, was acquired 
from hourly observations, and averaged over three-hour 
mm@eenvals., (Fig. 28) 

The boundary conditions, i.e., heat and momentum fluxes, 
epered were computed within the program. (Figs. 29, 30) As 
Will be discussed later, the average daily net heat flux 
forcing the Garwood and Mellor models over the thirteen days 
emer. 29 W/m=*2, (where positive is upward flux, i.e., the 
ocean was losing heat to the atmosphere) 

Additional information required by the Garwood model 
includes various coefficients; such as, the extinction 
coefficient and the drag coefficient. For the most part, they 
were set from previous model usage. The Garwood model con- 
tains two 'tuning parameters' (resulting in essentially two 
degrees of freedom) which had to be addressed prior to model 
Bemulations with each ‘new host point’ [Garwood, et.al., 
1982] The parameters, AM and AZR had already been set upon 
receipt of the model. (A tuning analysis is presently 
maaerway within the Ocean Prediction Laboratory. The procedure 
and results will probably be documented at some future time.) 

Sgeniiciemts, such as the diffusion coefficients, used 
in the Mellor level-2.5 model are ‘universal! constants, and 


SO required no alterations. 
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Mie jo bAristics 


pee ee LERMINATION OF THE OBSERVED VARIABLE MIXED LAYER DEPTH 

Oeeaaning mixed layer depths from DAS or FNOC digitized 
@rorites or even from the original XBT traces is no easy 
ask. Criteria employed vary greatly, and there is a ten- 
Bem@ey tO change one’s criterion with each new exception that 
Meosece Further complexity occurs with the desire to isolate 
new (shallow) mixed layer depths in addition to the numerous 
old (deep) ones. 

The new mixed layer is a layer formed over the old one 
due to the mixing being insufficient to reach the level of 
the old (previously formed) mixed layer. This may be the 
result of weak winds and surface heating. The old mixed 
layer will remain in place until turbulent diffusion of the 
interface or some other process dissipates it or strong mixing 
deepens the layer. The main reason for this is that the 
Second Law of Thermodynamics does not permit the unmixing of 
the fluid. An ‘excellent' example of an old mixed layer 
Meeorded during CCSI is depicted in Fig. 31. (This layer is 
also the new mixed layer for it has replaced the old one 
which would have existed at an earlier time.) An old mixed 
layer extending to approximately 63m with an overlying new 
mixed layer extending to approximately 41m is depicted in 


fe. 32. In Fig. 3, h, would be considered a new mixed layer 
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Syeriying the old winter mixed layer, hh. New mixed layers 
sem7Ooe depend solely on seasonal or synoptic changes. 
Numerous cases of diurnal mixed layers were noted wherein 
Bm@encooling of the surface at night caused the upper water 
Bemumn tO become unstable and consequently, mixed, forming a 
new mixed layer. Meanwhile, during the day, the surface is 
heated, possibly producing a warm shallow layer over the old 
mixed layer, 1.e., a new thermocline is formed. 

New mixed layers may also be formed under a large number 
Seeormer conditions. A storm, for example, will most likely 
Cause a new mixed layer to form which will quite possibly 
Be deeper than the old one. In this case, the new mixed 
layer supplants the old one, and becomes the old mixed layer, 
Ea . 

Eventhough the old mixed layers are of importance to Navy 
applications, for their role in surface acoustic ducting, the 
new mixed layers are the layers of importance in mixed layer 
modeling, i.e., that is the layer predicted and modified by 
some bulk (layered) models in use today. Thus, one must 
Med a suitable method for defining both of these significant 
features, as well as sea surface temperature, depicted in an 
MBI profile. 

Miestirst attempt to obtain mixed layer depths from the 
temperature depth profiles was made by computer comparison 
Seeehree sequential points on the trace for a 0.2C change. 


(The value 0.2C was chosen for that is the accuracy of the 
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Meee sts; more precisely, 0.12C-0.31C [Heinmiller, 1982]. 
It is the value presently in use in producing the TOPS/ 
TOPS-EOTS contoured mixed layer depth product.) This 
technique proved to be satisfactory except in cases of a 
surface thermocline, a very weak new mixed layer, no mixed 
layer, and an old mixed layer with a noticeable gradient 
(though still possessing the characteristic 'knee'). A 
Similar technique was used by Thompson (1975) with the 
additional technique of starting his sequential comparison at 
three meters depth to avoid the surface thermocline problem. 
With either procedure one would expect to obtain a depth, no 
Geweeer whether a mixed layer existed or not. (The depth for 
the 0.2C change obtained may be in the middle of a gradient.) 
such was not the case, however, due to the digitization 
mememe Used by FNOC wherein only depths of inflection points 
were listed. Comparison of the surface temperature with the 
memeperature of the first inflection point, the point at the 
base of the shallow surface thermocline, may yield a greater 
Mee 0 .2C change. This would in turn lead to 0.0m being 
@enoted the mixed layer depth when a very noticeable mixed 
layer might be apparent below the shallow thermocline or 
meecibly due to the response time of the XBT recorder. 

Mmigese possible deficiencies made it necessary to return 
to the original traces and logs. Each trace was examined and 
redefined in terms of temperatures and depths of the features 


of interest, for example, surface temperature, the temperature 
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memeaepth of the base of the surface thermocline, the 
memperpacture and depth of each sequential possible mixed 
layer, and the temperature and depth of a deeper point chosen 
Bomdertine the seasonal thermocline. Final mixed layer depth 
values, as well as other values, were chosen by comparison 

of the computerized analysis, the traces, and the XBT logs. 
The computerized analysis was given the most weight, followed 
by results obtained from inspection of the traces. The XBT 
logs were given the least emphasis due to their variation 
with the experience, knowledge, and judgment of the observer 
and the not always optimal conditions under which they were 


recorded. 


B. STATISTICAL CHARACTERIZATION 

The several variables derived from the XBT profiles by 
the techniques previously mentioned were used as a basis for 
a statistical characterization of each cruise. Basic pro- 
perties such as the mean, minimum, maximum, median, and 
standard deviation were examined. They are summarized in 
Table III. Two types of correlations coefficients were 
calculated, as well as Student-t confidence intervals, 
regression equations, and histograms and maps of SST, bucket 
temperature, old mixed layer depth (OMLD), old mixed layer 
depth temperature (OMLDT), new mixed layer depth (NMLD) and 
new mixed layer depth temperature (NMLDT). Only a limited 
selection of statistical quantities will be presented to 
convey the character of the upper ocean during the various 


eruises. 
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The first variables to be discussed are the SST's for 
the three cruises and the bucket temperatures from the first 
Meltse. A total of 105 XBT's were available from the CCSI 
grid for the statistical characterization. Both XBT SST and 
bucket temperature were recorded with a comparison revealing 
the XBT SST field was overall slightly warmer than the bucket 
Memeeature field, i.e., 0.4C. (Fig. 33, Table III) A 
regression equation or least squares linear equation of the 
Paeket temperature with respect to the XBT SST was calculated 
Memes: Y=.678 + .926(X1) where X1=XBT SST. The R squared 
value, a measure of how well the regression equation fits 
the data, was 79.6%. (100% is a perfect fit. The R squared 
value is defined by (100)(sum of squares due to regression/ 
total sum of squares.)) 

The two correlation coefficients computed were the Pearson 
Product moment coefficient and the autocorrelation function 
feeeerespect to distance. The former calculated for XBT SST's 
versus bucket temperatures, resulted in a large correlation 
coefficient, .892., though not ‘large enough! considering 
the fact that the values should be exactly the same (except 
for measurement noise). The spatial autocorrelation function 
was computed using binned XBT pairs. Out of 210 km, the 
radial resolution for XBT pairs within a bin was 195 km. 
Beyond 210 km, it was increased to 25 km. XBT SST for CCSI 
(Fig. 34) had a correlation distance of approximately 30 km, 
using a 0.2 correlation as a reference. The correlation 
distance for the bucket temperature was slightly less, l.e., 


27 km. (Fig. 35) 
100 
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During the latter two cruises, bucket temperatures were 
momeeeccorded routinely. The CCSII data set consisted of 97 
XBT's. As was to be expected, the SST values were much 
higher than those recorded for the March cruise. In fact, 
the minimum CCSII temperature, 12.9C, exceeded the maximum 
CCSI, 12.4C. The range of temperatures for CCSII, and thus 
the standard deviation in this case was quite large. A 
Memeo tor the significant range, 12.9 - 18.0C, was the 
presence of cold eddies and cool anomalies near the peak 
of summertime warming. 

mies autocorrelation functions for SST for CCSII were 
plotted with and without a linear trend removed, Fig. 36. 
This was done to get some idea of the existence and effects 
Sea crend in the data. Also shown is a 95% confidence 
region computed by +/- sqrt(2)/N where N is the number of 
Paris in the bin. (Using this N instead of the total N 
results in a more conservative confidence region.) Results 
indicate a slight increase of the correlation distance 
compared to the March cruise, i.e., 32 km for the detrended 
data and 38 km without the trend removed. 

The mean SST for CSSIII was still higher, 16.2C, as might 
De expected for later in August. The range, however, fit 
within that of the second cruise. This may be related to 
Mm@emeract that CCSIITI did not cover as large a region as CCSII. 
The correlation distance for this cruise was 40 km for the 


detrended data and 43 km without the trend removed. (Fig. 37) 
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This might indicate a more homogeneous area as compared to 
Biewrarger grid transited during CCSII. 

The new mixed layers for the March cruise ranged from 
0.0 m (when no new mixed layer was observed) to 82 m. (Table 
IfI) The largest Pearson products moment correlation obtained 
was for the NMLD's and OMLD's. This resulted due to a number 
Seeeases When the NMLDwas identical to the OMLD, i.e., in 
Cases of intense turbulence. 

For CCSII, 84 XBT's possessed a NMLD. The NMLD's were, 
however, quite small, approximately half of those noted in 
Meee trig. 38). CCSIII had 94 instances of NMLD's. The 
distribution in the histogram appears to be Raleigh in 
Meture. 

Histograms of the OMLD's for the cruises are shown in 
Fas, 39. Most noticeable is the shift toward shallower 
mixed layers with time. In fact, the mean depth for the March 
Cruise is a factor of two greater than the CCSIII mean OMLD. 
The Pearson product moment correlation coefficients for the 
OMLD in association with the majority of the other variables 
showed no dominant results. The largest coefficient obtained 
Was with respect to the NMLD; the probable reason for this 
was already discussed. For CCSI, the correlation distance 
was merely 15 km. (Fig. 40) For CCSII, both the data with 
and without the trend removed, have a correlation distance 
of only approximately 14 km. (Fig. 41) CCSIII possessed old 


mixed layers that were more highly correlated than those 
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previously examined. (Fig. 42) The correlation distance for 
mieedetrended data was approximately 35 km; without the 
trend removed the correlation distance was approximately 50 
km. The increase in the correlation distance from CCSII to 
CCSIII may be due to the fact that the former cruise covered 
a wider and more diversified domain than did the latter. 
Memeours of OMLD are included for CCSII and CCSIII to 
give an idea of the variability of mixed layers (a quasi- 
Eym@eptic depiction), and to bring to light a common problem 
Smcountered in trying to contour ‘real! data. The data must 
De placed on a regularly spaced grid prior to using contouring 
graphics packages. One would think the answer to the problem 
would be to merely interpolate the data to the new grid. 
Problems arise, however, from data set to data set depending 
on the noise or gradients in the data. Certain schemes seem 
to work well for some sets and not to work for others. The 
enclosed computer plotted contoured fields are the result of 
mieerpolation by trying to fit a five degree polynomial to the 
Mata. (rig. 43) The results do not have much worth as seen 
when compared with the hand contours. (Fig. 44) Only the 
locations of some of the minimum/maximum centers are retained. 
In the future, the technique of optimal interpolation should 
be applied to the data for improved interpolated values and 


Consequently, more appropriately smoothed contoured fields. 
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Figure 44, Hand contoured old mixed layer depths 
(meters) from CCSIII 
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Pee ooo lo Or OBSERVED AND SIMULATED 
Peiea. ATURE s TRUCTURE 

Mhe first run was of the Garwood model initialized with 
m@emlOors predicted profile for 1 August. The initial mixed 
layer depth was set at 20m, and the initial SST chosen was 
¥7.7C (the SST of the TOPS analysis profile). The 24 daily 
profiles were averaged so as to produce one 'representative' 
profile per day. The average profiles along with their 
Standard deviations compose Figs. 45-58. The standard 
deviation was nonzero only in the region at and/or above the 
base of the mixed layer, because the Garwood model does not 
alter regions below the furthest extent of the mixed layer 
base. For the 13-day cruise period, the Garwood model 
cooled the mixed layer by approximately two thirds of a 
degree, and deepened the mixed layer to approximately 37m. 
A substantial amount of entrainment occurred from 2 to 3 
August, with the simultaneous increase in wind speed (7.8 m/s), 
and decrease in the SST at 00Z. Thus, there was increased 
forced convective activity due to increased wind stirring and 
increased free convection due to surface ccoling, Significant 
Mixing also occurred over the time spans; 7 to 8, 9 to 10, 
12 to 13 August with some entrainment also occurring from 
Meco 12 and 10 to 11 August, but not as pronounced as that 
during the previously stated times. Not much occurred during 


the other time spans. 
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Figure 46. Average 1 August profile from Garwood 
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Figure 49. Average 4 August profile from Garwood 
Model. Zfnitialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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Figure 50. Average 5 August profile from Garwood 
model. Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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Figure 51. Average 6 August bDrofile from Garwood 
model. Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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Figure 52. Average 7 August profile from Garwood 
model, initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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Figure 53. Average 8 August profile from Garwood 
model. JInitialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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Pigure 54. Average 9 August profile from Garwood 
model Initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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Figure 55. Average 10 August profile from Garwood 
model. initialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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Figure 57. Average 12 August profile from Garwood 
model. initialized with FNOC analyzed 
profile (Lemperature in degrees Celsius) 
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Figure 58. Average 13 August profile from Garwood 
model. JInitialized with FNOC analyzed 
profile (Temperature in degrees Celsius) 
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The time-depth contours of temperature in Fig. 593 are 
another way of indicating the deepening and shoaling of the 
mixed layer depth. It is quite obvious that the depths 
below the mixed layer are not affected by the model. 

The TOPS profiles, on the otherhand, are affected at all 
levels (down to 500m). The profiles are more complex and 
harder to interpret, because of the variations occurring 
over the entire profile. The process of energy addition/ 
subtraction to the mixed layer is no longer as obvious as 
with a bulk model. Overall, changes ouee the 13 days, 
however, clearly show a substantial increase in the mixed 
layer temperature; the SST increased approximately one 
degree, while the mixed layer depth increased to approximately 
Mit rigs. 60, 61) (TOPS and TEOTS profiles for the other 
days may be found on the composite figures to be discussed 
in the last section.) 

A comparison of the contours for the TOPS predicted 
profiles (Fig. 62) with those of the Garwood model gives 
further indications of the difference in results of the 
Mr rusion and the 'slab' or integral models. For these 
contours to resemble those obtained from the Garwood model, 
the mixed layers would indeed have to be 'slab-like'. (The 
Somateaus’ located in the center of the contour field are due 
to repeated profiles, i.e., profiles for 6 and 7 August were 


met Obtained from FNOC.) 
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Figure 61. TOPS and TEOTS profiles for 13 August 
ge 2 (Temperature im degrees Celsius) 





seTtyodd SdOL JO pesodwos pleats Z/_ padnojuog 


SAVG 
O'*7t OLT O'et OIF OOT O6 O88 OL O9 OC OF OE Oe OF 


HUNOLNOD Z/AL DONG 





"Zg ednsty 


(W) Hida 


(cis 





Both of these sets of profiles and contours were compared 
with those obtained from daily averaged XBT profiles. The 
Variability of the upper water column as noted in the latter 
two sets of profiles was not as apparent in the series of 
averaged XBT profiles, as expected. The trend of the latter 
XBT's, averaged in space and time, had a tendency to wander 
Metre a bit. For example, the SST from 1 August to 2 August 
Jumped 1C with a corresponding increase of 0.3C at 400m. 
(Figs. 63, 64) Similar changes occurred for other profiles. 
meeecxample, the 13 August profile indicated a -1.3C change 
in SST and a -0.7C change in the temperature at 400m since 
Memrueust. (Figs. 65, 66) (Average XBT profiles for the 
other days may be found on the composite figures to be 
discussed in the last section.) Mixed layers and their tem- 
poral variations are not well represented by these profiles 
due to the spatial heterogeneity within the study domain. 

It takes a stretch of the imagination to identify a mixed 
layer in some of the profiles. Once identified, the mixed 
layers do not tend to exhibit the characteristics expected 
for mixed layer evolution, e.g., the deeper layers are not 
Meeessarily cooler. The characteristics seem to be more of 
eee random’ nature, i.e. they are affected by factors other 
mem those of mixing and heating. 

The use of average XBT profiles is somewhat of a contro- 
versial issue. One side holds that the averaging of XBT's 


has some worth in that it will cancel some of the spatial 
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Figure 63. 
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Figure 64. Average XBT and standard deviation plot 
for 2 August 1982. (Temperature in 
degrees Celsius) 
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-Pigure 65. Average XBT and standard deviation plot 
for 12August 1982. (Temperature in 
degrees Celsius) 
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end temporal effects. The other contends that the averaging 
of XBT's destroys the vertical structure, and, therefore, 
should not be the method by which one obtains a single daily 
profile with which to compare the model simulations. The 
latter believes that a median XBT (based on various cri- 
teria) would be a much better ‘measuring stick! for model 
Simulation success. Having gone through the abundant supply 
of XBT traces, comparing the averaged profiles with the 
series, the author is convinced of the latter viewpoint, 
though not contradicting that of the former. The corollary, 
then, is that the thermal structure is not robust over 


wy 


modest distances of a day's transit in a slow ship (ca. 3ms 
It is also suggested that greater care should be taken in 
grouping the XBT profiles in space-time for analysis. 
Indicated on the righthand panelsin Figs. 67 and 68 are a 
series of hypothetical XBT traces for two cases: one 
containing XBT's with pronounced jumps at the mixed layer 
bases, and the other containing XBT's with gradually sloping 
mixed layers. The lefthand panel of each figure depicts the 
average XBT trace for the series, as well as a median XBT. 
(The average profile is the shallower of the two.) From 
these it can be seen that averaging distorts the vertical 
structure, making the mixed layer depth shallower, and the 
temperature at the base of the mixed layer colder. 


fedian XBi profiles were almost as hard to define as 


were the mixed layer depths, especially over a substantial 
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memendand distance. XBT traces for a day varied from having 
fexcellent bulklike' mixed layers to having no mixed layers. 
Mifessot*s did not vary in any predictable pattern, i.e., 
sST's of XBT traces without mixed layers were in fact colder 
than the SST'’s of the traces with mixed layers. The transit 
of the ship across numerous fronts, jets, and eddies would 
be one explanation for these anomalies and ambiguities. 
Instead of obtaining a daily median profile, profiles 
for the same time each day were used as the ‘'observed' 
reality. The time chosen was 1800Z or 1000 local time to 
allow for the largest sample of XBT's, and to avoid periods 
Seemayor diurnal heating or cooling. Fig. 69 is the T/Z 
Goncour for the XBT's closest to 1800Z each day, approximately 
Within 30 minutes. There were no XBT's available for the 
designated time on 5, 6, 7, or 8 August. Therefore, the 
profiles for 5, 6 August were made identical to that of 4 
August. The profiles for 7, 8 August were made identical to 
the profile for 9 August. The 'wavelike' nature and the two 
vents present in the average XBI T/4Z contour field at 2, 12 
August are also located in the 18002 XBT T/Z contour field. 
Associated with these peaks or troughs, one might expect to 
find a change in the mixed layer due to fluctuations in the 
whole water column. Such in fact was the case for mixed 


layers of 2 and 1- August. 
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Returning to the results of the Garwood and TOPS model 
runs discussed thus far, one must keep in mind a number of 
important points. Recall that the Garwood model simulation 
shown was initialized with a profile from outside the study 
region and with inputs of varying degrees of accuracy. be 
was not known exactly how good the boundary conditions were, 
nor how sensitive that model would be to them. A sensitivity 
test was performed to decide whether the Garwood model was 
deepening the layer due to the 19.5m FNOC winds being too 
strong. This was done by decreasing the drag coefficient 
by 10%, equivalent to decreasing the FNOC winds to 95% of 
their original magnitude. As can be seen from profiles for 
mead 13 August 1982 (Figs. 70, 71), the effects are not very 
substantial. 

In addition, it should be remembered that the TOPS model 
was initialized and run at a point outside of the cruise 
grid. It is not known what values TOPS received as initial 
or boundary conditions, only that it received its atmospheric 
forcing from the atmospheric primitive equation model and 
NOGAPS, and that eight XBT's were used in the TEOTS analysis. 
Finally a slightly warmer SST is expected for the FNOC 
Ppadpoint, because it is further offshore than the study 
domain. According to the previous fortnightly mean, 
however, the difference should only be about 1C. 

To make simulations more consistent within the domain, a 
second run was performed using XBT 96 to initialize the 


Garwood model. This XBT, deployed on 4 August, was chosen 
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Figure 70. Daily average and average tuned Garwood 
profiles: 1 August (Temperature in 
degrees Celsius) 
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Figure 71. Daily average and average tuned Garwood 
profiles: 13 August. (Temperature in 
degrees Celsius) 
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because its location was resampled on 10 August with the 
mien Of XBT 216. (Fig. 72) (It was hoped that by choosing 
XBT 96 to initialize the Garwood model, its mixed layer 
evolution could be compared and evaluated with the data 
Seeained from XBT 216. The initial SST used was that of 
XBT 96. An initial mixed layer depth of 20m was "assumed". ) 

In addition, to avoid comparing simulations of two 
differently forced models, the Mellor level-2.5 model was 
obtained and modified to produce simulations using forcing 
identical to that used in the Garwood model. Profiles 
Pesulting from the model runs (Figs. 73-85) along with the 
applicable 1800Z XBT profiles, TOPS predicted profiles, and 
daily averaged profiles are examined. (Notes: A) repeated 
TOPS, 1800Z, and average XBT profiles during in-port period; 
B) TOPS profiles correspond to 1600 local time, il.e., 
afternoon heating exists; C) Mellor profiles are not daily 
averages.) 

Based on the overall change in the profiles, Garwood and 
Mellor profiles cooled and deepened the mixed layer. The 
profiles seemed to be quite consistent with the observed 
data by 4 August and even more so by 10 August. The cooling 
trend in the XBT's seemed to be consistent with that apparent 
from XBT 96 to XBT 216. However, caution is required, 
because the role of the horizontal and vertical advection 


Mas not been assessed. 
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Figure 72. XBT's 96 and 216 (Temperature in 
degrees Celsius) Launched at approximately 
the same location on 4 and 10 August, 
respectively 
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Figure 73, August 1 composite of simulated and 


observed profiles 
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Figure 74, August 2 composite of Simulated and 
observed profiles 
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Pigure 75. August 3 composite of simulated and 
observed profiles 
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Figure 76, August 4 composite of simulated and 
observed profiles 
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Figure 78. August 6 composite of simulated and 
observed profiles. 
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Figure 79. August 7 composite of simulated and 
observed profiles. 
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Figure 80. August 8 composite of simulated and 
observed profiles 
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Figure 81. August 9 composite of simulated and 
observed profiles 
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Figure 82. 
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August 11 composite of simulated and 
Observed ‘profiles. 
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Figure 84, August 12 composite Of simulated and 
Observed profiles. 
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Figure 85. August 13 composite of Simulated and 
observed profiles. 





Mites lYPS profiles, though not expected to resolve the 
small scale veriability of the boundary current region, on 
the otherhand, started out warmer than the other profiles and 
indicated additional warming with time. The fact that the 
gridpoint used was further offshore than the cruise region 
caused one to expect a slightly greater surface temperature 
memeene FNOC profiles. That the forcing averaged 184 Watts/ 
m**2 per day was, however, not anticipated. This heating 
is, however, believed to be part of a larger problem, viz., 
a heat bias in the atmospheric PE model. Heats biases of 
varying signs have been previously identified by Elsberry, 
Seem (1979), Budd (1980), and Warn-Varnas, et.al. (1982) 
and shown to possess temporal and spatial scales of a month 
and tens of thousands of kilometers, respectively. Warn- 
Varnas, et.al. (1982), for example, found that a limited- 
mg@eeeversion of TOPS initialized and forced with FNOC fields 
revealed considerable skill during 60-day runs in comparison 
to persistence if a bias of about 150 ly/day was removed 
from the daily net surface heat flux of the domain. This 
correction proved to yield consistent results with the heat 
fluxes computed directly using atmospheric observations 
acquired by ships-of-opportunity. (Clancy, et.al., 1980] 

In another study, comparison of the FNOC heat fluxes with 
those calculated by Nate Clark indicated a cold bias believed 
to be mainly the result of excessive latent heat loss. The 


problem of heat bias is extremely evident in such data sparse 
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Seeds as the high latitudes. POrecxample, I hOls Sol's in 
Biemwoea Of Okhotsk, the Bering Sea, and the Labrador Sea 
pave in the past exhibited a warm bias, running as much as 
6C above those of conventional EOTS (essentially 
representative of climatology within these areas). [Clancy, 
et.al, 1982] 

In the simulations presented, the initial conditions and 
forcing for the first two days were supplied to TOPS by the 
PE model. On 3 August, NOGAPS replaced the PE model. 
However, any spurious heat available in the PE model entered 
mero tie simulation prior to bringing NOGAPS online. Once 
introduced into the run, it is unlikely that the heat could 
be removed without much mixing or surface cooling (due to 
the Second Law of Thermodynamics). NOGAPS' contribution to 
the spurious heating is unknown, though it's clear that the 
heating trend was not corrected over the short run time. 

(It is as yet unknown whether NOGAPS has a heat bias of any 
kind.) 

Examining the TEOTS SST anomalies (TEOTS - climatology) 
for some of the cruise days as well as for a number of days 
following the cruise, and following the ‘return to climatology' 
adjustment installed in TOPS on 20 September (Table IV) 


revealed an excess of 2-3C in the study domain as compared 


lmorconventional EOTS. 
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TEMPERATURE ANOMALIES 


a_i 


DATE TEMPERATURE ANOMALY (C) 
i es se C(‘(C‘Ct;tS 

1 SIS EOrS 

820730 2 

820803 NOGAPS Online 

820804 Z = 

820808 Z 0 

820811 3 

820814 2 

eZzOorg . 

S206 27 Z 

820828 system Restart 

e040 Hive 

20920 1.4 

eZ0g23 Ess, 


*Return to climatology adjustment installed 


The warmer SST's along with the an identified tendency 
for the turbulent parameter scheme used in TOPS to underpredict 
the MLD's [Martin, 1982] might be important in explaining the 
trend in the TOPS' profiles. The difference between various 
model profiles with depth is believed to be the result of the 


strong connection TEOTS has with climatology. 
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To make a strict comparison of model results (assuming 
ideally that the model forcing is the same, as it should be), 
one should interpolate between FNOC gridpoints to the 
Meeations of XBT 96. 

It is felt , however, that the heating difference is not 
fotally a function of distance. Erroneous heat fluxes, 
strong dependence upon climatology, and advection (associated 
with mesoscale variability) are thought to play important 


mMebes in the results. 
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V. SUMMARY AND CONCLUSIONS 


The Garwood and Mellor level~-2.5 simulations resulted in 
cooling and deepening of the mixed layer over time, and 
exhibited a likeness to the observations. The TOPS model, 
receiving different forcing, produced warmer profiles with 
a substantial warming trend over the 13 days. 

The exhibited differences of the TOPS profiles with 
respect to the profiles produced by the other models are 
believed to be more than a result of location. They are 
believed to be due to a) a heat bias in FNOC's atmospheric 
PE model, b) too much dependence upon climatology, and 
Syeadvection and the subgrid scale variability existing in 
the region. 

miemautocorrelation functions for SST and MLD indicated 
a Maximum correlation distance of 35 km. However, TOPS 
uses a standard FNOC grid with spacing 381 km at 60N. 

Hence, though interpolation should be carried out when looking 
at regions between the gridpoints, it seems highly unlikely 
to be able to recover the small scale variability dominating 


the region. 
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VI. RECOMMENDATIONS 


First of all, numerous recommendations are in order 
concerning the data acquisition on future cruises. Most 
importantly, all easily observed data should be recorded 
regardiess of the immediate mission of the cruise. (e.g., 
wet bulb temperature) The equipment should be in a working 
Status, and be sufficiently backed-up with manual logs. 
Underway analysis should be carried out in a consistent 
manner wherever possible to avoid excessive reworking of the 
Semeaeashore. A thorough scientific crew indoctrination 
period is, therefore, in order. 

From the modeling aspect, many more evaluations of TOPS 
are necessary under numerous varied conditions prior to its 
full acceptance by the scientific community. It is important 
that the TOPS products be interpolated to the observed 
position and time of the study domain prior to any further 
evaluations of the model's performance. Design tests should 
be carried out with control over the variables. Independent, 
direct measures of atmospheric forcing should be made to 
drive the research and operational models. (Simulations 
should be produced using both sets of forcing values). 
Additionally, it would be of interest to attempt to predict 
SST and MLD from one realization to the next, and to examine 


the changes in the profiles in the (2,x,t) domain instead of 
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merely the (z,day) domain. All of these tests and comparisons 
should also be performed on the 60 km and 20 km grids of the 
future. 

In addition, other operational models and their products 
mpe becoming more important both in the Navy and in the 
Civilian community. With their growing availability, so too 
must grow the amount and quality of verification testing 


accorded them, for they are an important asset of the future. 
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APPENDIX A 


As DATA PROCESSING PROCEDURES 

Working with a new data set is rarely an easy task. The 
ieaemust be edited and processed into a usable format. 
Here, the editing procedures used to produce the five 
required inputs for the version of the Garwood model at hand 
e@escescribed in more detail to better assess the model 
inputs and results. 

The procedure applied to the DAS data was as follows. 
First, the data on magnetic tape was dumped onto mass storage 
and filtered for blocks of zeros and gross errors in lati- 
tude and longitude. Following this, the data from the 
entrance and exit tracks to the cruise grid were deleted so 
aomeeo cCOnfine the data to the study domain. Finally, the 
finer errors in temperature, wind speed, etc. were deleted. 
The resulting air temperature data were then averaged over 
three-hour intervals to yield the required eight daily air 
temperatures. Averages were computed over bins containing 
from 6 to 299 values with the maximum standard eye hoe I Ope. 
any resultant value being 1.67C and the minimum being 
Mmaedesc. Only three of the averaged values possessed 
standard deviations over 1C. Procedures used to obtain 
Meesamnes Values within cruise days were based on the 
assumptions described previously. For days in port, the 
air temperatures were held at a representative constant value. 


(This could be improved upon with reasonable means to predict 


gaps.) 
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The dewpoint temperatures were erroneous due to equipment 
failure. The required eight values of dewpoint per day were 
calculated by assuming a relative humidity of 90% and using 
Gim-mprocedure described in Bolton (1980). This value for 
relative humidity was chosen to correspond with the 
intermittently hazy/foggy conditions. 

Cloud cover indications from the ship's log were either 
in tenth of sky cover or in descriptive terms. The first 
eeeeeaken was to convert the descriptors to tenths of sky 
Fever using Bowditch's American Practical Navigator and the 
experience of the R/V ACANIA's first mate. The tenths of 
sky cover were then converted to octals,and averaged 
over three hour intervals. Missing values during cruise 
days were not a problem as they had been with the air and 
dewpoint temperatures. The cloud cover during the in-port 
days was held at a constant five tenths. 

Winds could have been obtained in two ways. The DAS 
recorded winds could have been converted to true winds and 
averaged as was done with the air temperatures. The other 
way winds could have been obtained was by interpolating the 
daily analyzed TEOTS 39.05N, 127./75W wind to the necessary 
three-hour interval. The TEOTS analyzed winds over both 
@e>'t and CCSIII are depicted in Figs, 25 and 26 with the 
Meter time series containing the interpolated eight three- 
hourly values. (The symbol denotes the value at 00 Zulu) 


mime latter method was chosen for convenience as well as to 
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Facilitate the comparison of the results of the TOPS, Mellor 
Mewes =2.5, and Garwood models. Some of the 'CTD station 
winds’ (already assumed to be true winds due to the ‘zero! 
speed of advance of the ship) were plotted against the FNOC 
archived winds to observe whether or not the latter could 
‘reasonably' be applied to the cruise area. The CTD winds 
were found mainly to fluctuate around the FNOC winds except 
for the first day and a half when the 'CTD winds! were 
Sreeeeaed by the FNOC winds. The FNOC winds were (Fig. 86) 
employed in subsequent runs. (The triangles represent 'CTD 
winds’ recorded along with suspect data; leading to a 
possible suspicion in the winds.) A sensitivity analysis was 
done to assess the importance of the wind stress as a bound- 
ary condition, and in so doing, ascertain the magnitude of 
possible problems introduced by using winds recorded at 
19.5m vice 10.0m. The procedure and results were described 
iieehapter IV. 

The final input required was SST. The SST's were avail- 
pees trom a number of sources. DAS recorded boom temperatures 
could have been averaged. This route was not chosen, because 
of the possible problems posed by the boom varying in its 
depth in the water (interaction with the ship's wake). 
Alternatives included the use of the bucket or XBT surface 
temperatures, averaged or singly, at the required times. 

For the initial runs, the XBT SST's recorded by DAS 


were used: this was believed to be the least subjective 
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approach. The bucket temperatures were subject to the 
experience of the watchstander, as were the XBT SST's 
recorded in the XBT logs. The validity of the latter were 
found to strongly depend upon the existence of a surface 
Biemmectine, either 'real' or possibly due to the initially 
nonlinear slow response rate of the XBT recorder. The XBT 
imemeeratures were taken from those launched closest to the 
required input time. (usually within 30 minutes--sometimes 
linear interpolations were performed.) Averages were not 
made for the three-hour transit of the ship across possible 
meOnts . 

Finally it was necessary to input an initial temperature 
profile into Garwood and Mellor level-2.5 models. Temperatures 
were interpolated to intervals of one meter to a depth of 
200M. The first runs were made using an interpolated TEOTS 
profile corresponding to 1 August. (The original TEOTS 
analyzed temperature profiles consisted of 13 temperature 
depth values. Surface values were not available; the tempera- 
ture at one meter was used at the surface.) Subsequent runs 
were made using XBT 936 as the initial temperature profile for 
the Garwood and Mellor level-2.9 models. This XBT was chosen, 
because of its proximity to the FNOC gridpoint of interest. 
More importantly, however, 1t was chosen because of its common 
Meet iOm with an XBT, 216, launched six days later. XBT's 96 
and 216 were launched on 4 and 10 August, respectively, thus 


Spanning a part of both cruises. 
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APPENDIX B 


TABLE V 
eee SOUL TON OF FORS 


VERTICAL GRID SPACING USED IN TOPS 


Level Depth (m) 
i ass 
2 se 
3 LASS 
si ae 
5 Laer 
6 SAS 
J yO 
8 S020 
g 62.0 

10 ivomeee) 
Wea a3) 0) Be 
ie? enone) 
iS 15080 
14 PAONONe 0 
19 S000 
16 aekemae 
ay SUG Go 


TOPS is run on the Cyber 205 computer on a standard 
63x63 FNOC global grid. Values of T, S, U, and V are defined 
with this horizontal resolution at each of the seventeen 
levels listed above. All spatial derivatives are centered 
in space. The horizontal advection components are defined 
on a staggered grid; they are forward-differenced in time. 
Vertical eddy fluxes and W are defined midway between these 
levels. The former are differenced backward in time, while 
the latter (as well as the Coriolis term) is differenced 


trapezoidally in time. 
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Additionally, TOPS is quasi-three-dimensional model. 
(though a nonadvective version is kept on the Cyber 175 as a 
backup.) Advection, as previously mentioned, is by 
climatologically averaged geostrophic currents and 
instantaneous wind drift. The geostrophic currents are 
updated monthly using density fields computed from monthly 
temperature and salinity climatology and the thermal wind 
Equations integrated upward from a predetermined "level 
of no motion", The geostrophic currents are nondivergent. 
The sole contributor to the vertical advection component 
is, therefore, the instantaneous wind drift via the Ekman 


eweeion mechanism. 
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